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Abstract 

The phase diagram of the Gross-Neveu (G-N) model in 2+1 dimensions as a function of chemical 
potential and temperature has a simple curve separating the broken symmetry and unbroken 
symmetry phases, with chiral symmetry being restored both at high temperature and high density. 
We study, in leading order in the 1/N expansion, the dynamics of the chiral phase transition 
for an expanding plasma of quarks in the Gross-Neveu model in 2+1 dimensions assuming boost 
invariant kinematics. We compare the time evolution of the order parameter (mass of the fermion) 
for evolutions starting in the unbroken and broken phases. The proper time evolution of the order 
parameter resembles previous results in the 1+1 dimensional G-N model in the same approximation. 
The time needed to traverse the transition is insensitive to \i. 
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The phase structure of QCD at non-zero temperature and baryon density is important 
for the physics of neutron stars and relativistic heavy ion collisions. The Gross-Neveu 
model [|IJ has been a quite successful toy model for understanding the general features 
of chiral symmetry breakdown. In a previous paper B we studied the dynamics of the 
chiral phase transition in leading order in 1/N in the 1+1 dimensional Gross-Neveu model. 
Although the phase structure in that approximation mimics that of two flavor QCD, in 1+1 
dimension the exact theory only has chiral symmetry breaking at zero temperature. In 2+1 
dimensions, however, the Gross Neveu model with Z 2 chiral symmetry breaking does indeed 
have interesting phase structure in the exact theory as shown by lattice simulations |§ . The 
lattice simulations also show that most of the phase structure can be understood from the 
1/N result. The phase diagram obtained in leading order in 1/N is also qualitatively correct. 

The Gross-Neveu model in 2 + 1 dimensions was studied in the large-N expansion exten- 
sively in Refs. (0-0). The Lagrangian for this model is given by 

C = -t%-fd^i - ^ 2 (^) 2 , (1) 

where the index i = 1 . . . N. Here the 7 M are the conventional four dimensional Dirac gamma 
matrices, but with \i = 0, 1, 2. By introducing the auxiliary field a = —ig 2 ^>S>i/N one can 
rewrite this Lagrangian in the alternative form 

Na 2 

C a = -M i [fd ll + <T\9 i - irT . (2) 
T = S eff = J d 3 xC 2 [V d ,a d } -iNtr In S-^aa}. (3) 

where 

S- 1 [a] = -i[^d, + a}5{x-y). (4) 

Letting = N then scales out of the action. The effective potential at finite tem- 
perature and chemical potential for this model has been discussed in detail by Rosenstein, 
Warr, and Park and to leading order in large N the effective potential is expressible as 



-2 



d 2 k 



u k + - ln(l + exp[-(3(uj k - /i)]) + ln(l + exp[-/3(u; fc + //)]) 



(5) 



where uo\ = \k\ 2 + a 2 . The derivative can be evaluated explicitly in 2+1 dimensions 

dV_ !L-2a ( d2k 

da g 2 J cj fc (27r) 2 

a J u k (2n) 2 \l + exp[-p(u k - /i)} 1 + exp[-/3{u k + /i)} J 
The temperature dependent part is finite and is given by 

4- ln(l + 2e~ /3W cosh^/3 + e' 2 ^) 

Tip 

The vacuum contribution has a linear divergence, which is related to coupling constant 
renormalization. The vacuum structure in the cutoff theory can have two phases depending 
on whether the bare coupling is less than or greater than the critical coupling: 

_ 2 r^E fp^^_ Mm 2Ae = Am m 

9mt i-A B (27r)3p| J-A M (2n) 2 2\p\~ tt 2 U 

Since we are only interested in the phase where there is broken symmetry in the vacuum 
sector, we can use the vacuum gap equation to renormalize the leading order in large- 
N effective potential. The value of a at the minimum of the potential for the vacuum 
determines the Fermion mass 

^\a=m f = = ^ + trS^J \ a=mf -y m f = -g 2 trS[m f }. (8) 

This leads to the well known gap equation 

1_ M 2 p 2 dp 1 = A-\m f \ 

g 2 Jo (2ir) 2 p 2 + m 2 f tt ' 1 ' 

where we are using a Euclidean cutoff. Using the gap equation, we therefore have that 

and we obtain in the vacuum sector 



dv = _ 2(y r dp E p 2 E 

da J 7r 2 

Integrating we obtain 



p\ + a 2 p 2 E + m) 



\\a\-\m f \). (11) 



TT 



v= lM_rn£\ (12) 



Tt\ 3 2 

Restoring the effects of finite temperature and chemical potential we have 
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FIG. 1: Plot of the line of phase transition given by Eq.(17). The y-axis is the chemical potential, 
and the x-axis is the temperature. In this plot m/ = 1 and fee = 1- 
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1 / 1 cr | 3 m/<T 2 " 
tt l~3 2~ 



2 f d 2 k /, p e+„ 

J J I [ + e_ ~ ] + [ + e ~~ )] 



and 



^IL = a J\ a \_\ m \ + \ ln[ i + 2e -^l cosh[///3] + e - 2 ^l^l] 
acr 7T V p 



The finite temperature effective mass is the solution of 

o- = rn/ - - In [l + 2e~ /3H cosh[/i/3] + e" 2/3|<T| 
P 



(13) 
(14) 

(15) 



If we approach the line of phase transitions from the region in which chiral symmetry is 
still intact, we can require 



d 2 V f 



eff i 



da 2 



CT=0 



0. 



(16) 



which straightforwardly leads to the relation that defines the line of phase transitions, 

1 1 n e m ^ 

m f = - ln[2(l + cosh[/i/3])] - fi = - cosh^— - 1] (17) 

Using this result the phase diagram of this model is easily constructed and shown in Fig. ([[]). 
The phase transition described by V e ff being second order. The exact phase diagram is found 
in and is quite similar. 



DYNAMICS OF THE PHASE TRANSITION 



In order to investigate the dynamics of the phase transition in a regime relevant to Rela- 
tivistic Heavy Ion Colliders, it is convenient to switch to curved space-time coordinates and 



assume boost invariance for the longitudinal expansion of an initially highly Lorentz con- 
tracted plasma of quarks. This type of assumption has been verified recently in the Hartree 
approximation in scalar field theory (see pi). Thus we will assume that the expansion 
possesses hydrodynamical scaling x — vt |TD[ , where x is the collisional axis. 



The change of coordinates to a boost invariant frame is given by x = rsinh?7 and t = 
t cosh f]. The corresponding Minkowski line element is 

ds 2 = ~dr 2 + r 2 dr] 2 + dy 2 = (-1, r 2 , 1). (18) 

Vierbeins are defined by = V^V^rjab, where T) a b = (—1, 1, 1) is the Minkowski metric. An 
obvious choice for the vierbein is 

v; = (i,r,i), V2» = (i ) i,i)-.dety = v=^ = r. (19) 

The action for our model in general curvilinear coordinates is 

2 

STM = Jd 3 x dety^^V^ + ^V^)^-^^-^), (20) 



where 7 M = 7 a V^, V M = <9 M + T^, and is the spin connection given by (11], [12 



= VK,(W + r; x V b x ) and Z ab = i[ 7 a , 7 b ]. ( 21 ; 



The Christoffel symbols have the usual definition, = (d\g^ + d^g\ — d u g^)/2. The only 
nonzero spin connection is T v = — |7°7 1 . From Eq. (|20|) , we find that the equation of 
motion for \1/ is 

(7 M V„ + a)* = 0. (22) 
Writing this out more explicitly, we obtain 

( 7 2 <9 y + 7 °(<9 r + + ^d v + <x)V = 0. (23) 
Zr r 

If we rescale the equation by letting \1/ = &/y/r and r = -^e u , we obtain, 

(rj 2 d y + ^du + 7 1 <9„ + a)$ = (ffd„ + cF)$ = 0, (24) 

where a = err. 

We expand the Fermion field in terms of Fourier modes 

$>) = j Mj&^fa)*^ (25) 



where k ■ x = k v r] + k y y. Note that the dimensions of b and d are in fact y/r, and therefore, 
it follows that $ ± is dimensionless. These fields obey the canonical quantization conditions 

^ a (r,x),^(r,y)] + = 5 a ^ 2 \x-y), (26) 

where a, (3 correspond to the component index and run from 1 to 4. We take the standard 
anti-commutation relations for the creation and annihilation operators: 

[b s (k),bl,(q)} + = [d s (k),dl,(q)} + = (2Tt) 2 5 ss ^ 2 \k - q). (27) 

The orthonormality condition for each k is 

(*iU«)) t *iU«) = ^-U u ))^-kM = (28) 

The equations of motion for the mode functions are 

(t^rky + 7 °<9 U + + a) = 0. (29) 

This is the equation that we will evolve in the conformal time (u) for our simulations. 
Explicitly, the equations in component form for the each spin are 

(d u - icr)$M + ( k v ~ iTk y )<S> 2tS = 
-(d u + ia)$ 2 ,8 + (K + iTk y )<$> 1)S = 
-(d u + ia)®3,s - (k v - iTky)§4,s = 

(d u - icr)$ 4iS - (kr, + irky)^ = 

(30) 

where 1 — > 4 denote the component. 

For the initial conditions, we take a complete orthonormal set of solutions of the Dirac 
equation for a fixed initial value of the mass equal to cr(0), 

K s {u Q )=yl s = u ks e-^ 



with 



U ks 



^Au k (u k + cr(0)) 



(-i^k + m) Xt (32) 



1 



V-ks 



i^K + ^(o) ) x. 



(33) 



'4u! k (u k + cr(0)) 

where A;' 1 = {tDfc(O), k v , rk y }, oj\ = k 2 + cr 2 . In what follows we use the explicit representation 
given in [II, 



and 



12D = ?(T 3 



i ' 
.0 -i. 



a u 
lAD = 12D 



il2D = °"1 



I N 
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/o r 

U o, 



>7 2 D = °"2 
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1 o . 



-1,1,1,-1). 



We choose the xV 1 as the positive and negative eigenvectors of 

f°\ ( \ ( l \ / 1 \ 
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The orthogonality condition for these eigenvectors is 

( X 0l )tx°j = 28 i5 8 ss , 

where i = {+,—}• We also have 

^l xf = ±X° S ± , = xg, and *7 2 X?S = ±xS- 



(34) 

(35) 
(36) 



(37) 



(38) 



(39) 



Effective Mass 

In lowest order in large N, a is the effective mass of the Fermion. Its value determines 
the phase one is in, and it is the order parameter for this model in this approximation. Let 
us now obtain a renormalized equation for a. If we symmetrize a 

ff = ^<[*t >7 °«r]>, (40) 

and use \I/ = 3>/\/t and Eq. (|25|), we obtain after using the canonical anticommutation 
relations 

a = 

Y T J ^E^ 1 " 2JV + (fc))($it.(t*))V$+.(«) + (2J\L(*) - l)($ M (n))V$ M ( M) ). 

(41) 
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So the equations we need to solve are Eq. fl3~0|) for ($^ s (w)) and Eq. (|41| ) for cr. To simplify 
the update we realize that since 

(*t(«))W^>) = -( $ fc »)V^-, s («) (42) 

it is only necessary to solve the Dirac equation for the positive frequency modes. As it stands 
this equation has a linear divergence. The most convenient way of handling the divergence 
is to use the definition of the unrenormalized coupling constant (suitably written in terms 
of mode sums) to relate the bare coupling to the momentum cutoff and the renormalized 
mass. That is one uses the lattice version of 

- = 2 f" ^ , 1 (43) 
9 2 J— A (2tt) 2 Jp2 + m f 1 ; 



to define the value of g used in Eq. (f41|) . 



At the initial time Uq, Eq. (41) reduces to the equation for the self-consistent mass at the 
initial temperature and chemical potential. Using the intial time adiabatic wave functions 
for /, we find for each spin 



Whence, 



(<^ S M) V$t>o) = ±-- (44) 

, ._2g 2 o(u )rdk y dkJl-N + (k)-N^k)\ 
t J (2vr) 2 \ u k ) 



where we should recall that u>k = yk ■ k + cr 2 . 

Using the vacuum equation for the gap equation evaluated at a = rrif one obtains the 
finite equation for the initial conditions: 

2tt f dk y dk v ( N+(k) + N.(k) \ 
<7(u ) — rat H / / ' , , 46 

k 1 ^ \^k-k + o{u y) 1 ; 

where we have used rk z = k v . In dimensionless form this is written as 

. _ n r d 2 k (N + (k) + N4k)\ 

where rhf = rnif and d 2 k = rdk y dk v . This is the same equation as the equation we obtained 
from the finite temperature effective potential Eq. fll5|). 

The phase transition is quite similar to what we found in 1+1 dimension. We have chosen 
three points above the line of phase transitions to study the effect of chemical potential on 
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FIG. 2: Plot of a vs. u for the second order phase transition. In this plot the dashed line is for 
the initial conditions To = 0.3 and /io = 1.0. The solid line is for To = 0.6 and no = 0.8 We chose 
an initial value <j{uq) = l.e — 5. 

the proper time evolution of the transition. We used \i = 1.0, T = 0.01, \i = 1.0, T = 0.3, 
and jj, = 0.8, T = 0.6. In order to break the symmetry it is necessary to give a a small 
initial value. For smaller values than a(uo) = l.e — 5, there was no change in the onset or 
steepness of the transition. However for zero or smaller (t{uq) numerical errors could lead to 
the final value of a being ±1. We found that for each of these initial conditions, the phase 
transition occured at almost exactly the same value of u and progressed similarly after that. 
A comparison of two of these simulations using \x = 1.0, T = 0.3, and fi = 0.8, T = 0.6 is 
shown in Fig. (0). The result for /i = 1.0, T = 0.01 is quite similar. 

If we instead started in the broken symmetry phase and studied the time evolution we 
would obtain the behavior (for \i — .6, T — 0.5 ) shown in Fig. (H). We see now that for 
this case the proper time needed for the system to relax to the vacuum is shorter, and the 
transition is not sharp. 

FC and VS thank the Santa Fe Institute for their hospitality. VS thanks the National 
Science Foundation for financial support. We also thank Wenjin Mao for computational 
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FIG. 3: Plot of a vs. u for no phase transition. In this plot To = 0.5 and jiQ = 0.6. We chose the 
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